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ABSTRACT 

We have analyzed a uniform sample of 16 evolved H II regions located in a 2°x2° Galactic field centered 
at {l,b) = (30°, 0°) and observed as part of the Herschel Hi-GAL survey. The evolutionary stage of these 
H II regions was established using ancillary radio continuum data. By combining Hi-GAL PACS (70 /itm, 
160 yum) and SPIRE (250 /xm, 350 /xm and 500 /xm) measurements with MIPSGAL 24 ^m data, we built 
Spectral Energy Distributions (SEDs) of the sources and showed that a 2-component grey-body model is a 
good representation of the data. In particular, wavelengths > 70 /im appear to trace a cold dust component, 
for which we estimated an equilibrium temperature of the Big Grains (BGs) in the range 20 - 30 K, while for 
A < 70 /Ltm, the data indicated the presence of a warm dust component at temperatures of the order of 50 - 90 
K. This analysis also revealed that dust is present in the interior of H H regions, although Ukely not in a large 
amount. In addition, the data appear to corroborate the hypothesis that the main mechanism responsible for the 
(partial) depletion of dust in H II regions is radiation-pressure-driven drift. In this framework, we speculated 
that the 24 /im emission which spatially correlates with ionized gas might be associated with either Very Small 
Grain (VSG) or BG replenishment, as recently proposed for the case of Wind-Blown Bubbles (WBB). Finally, 
we found that evolved H II regions are characterized by distinctive far-IR and sub-mm colors, which can be 
used as diagnostics for their identification in unresolved Galactic and extragalactic regions. 
Subject headings: Galaxy - H II regions - dust 



1. INTRODUCTION 

H II regions are sites of ionized gas surrounding OB asso- 
ciations and represent powerful laboratories for understand- 
ing the mechanisms regulating massive star formation. Their 
evolutionary sequence spans different phases, from the hyper- 
compact (HCH II) phase, in which the newly-formed star 
is still fully embedded in its dust cocoon, to the evolved 
one, when most of the natal material has been ejected from 
the birth site. Intermediate phases are represented by ultra- 
compact (UCH II) and compact H n regions. Each of these 
stages is characterized by critical values of electron density, 
n^, and linear diameter (Kurtz 2005): HCH II have extremely 
high electron density (n^ > 10*' cm"-') and small sizes (diame- 
ters < 0.05 pc), UCH II and compact H II regions have lower 
Tie (Ug ~ 10^ - 10'* cm~^) and larger diameters (< 0.1 - 0.5 pc), 
while evolved H 11 regions present even lower electron den- 
sities (Ug ~ 10^ cm"-') and diameters of the order of several 
pc. 

The emission properties also vary with the evolutionary 
stage. The radio emission is dominated, in all cases, by free- 
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free emission associated with ionized gas. However, HCH II, 
UCH n and, occasionally (depending on n^) compact H II re- 
gions, have an optically thick spectrum for frequencies from 
a few to several GHz. On the contrary, the radio spectrum 
of an evolved H II region is optically thin already at ^ 1 
GHz. Both young and old sources present radio recombina- 
tion lines (RRL), although HCH II are characterized by sub- 
stantially broader lines, likely due to pressure broadening and 
gas bulk motions. All evolutionary phases, with the exception 
of HCH II (Kurtz 2005), are very bright at FIR and sub-mm 
wavelengths. This emission has often been ascribed to the 
presence of dust grains. However, it is not clear whether dust 
is directly associated with H II regions or rather with their sur- 
rounding Photodissociation Regions (PDRs, see HoUenbach 
& Tielens (1997) for a complete review), since the histori- 
cal lack of resolving power in the existing IR data has long 
prevented us from distinguishing these two components. This 
elusive separation has led many studies to use large photomet- 
ric apertures which often included both the actual H II region 
and the PDR, with the consequence that the derived proper- 
ties were a global average over intrinsically widely different 
regions. In fact, while an H II region is by definition the region 
in space where the gas is fully ionized, the PDR is mainly neu- 
tral, characterized by an A^ > 1. In addition, the intensity of 
the radiation field, as density, is not the same inside or outside 
the H II region, as both scale with distance from the ionizing 
OB association. 

This situation has dramatically changed with the release of 
high-sensitivity, high-resolution radio (MAGPIS, Helfand et 
al. 2006) and IR (GLIMPSE, Benjamin et al. 2003; MIPS- 
GAL survey, Carey et al. 2009; Hi-GAL, MoHnari et al. 
2010) data, which allows the clear identification of an H n 
region from its PDR, thus offering the unprecedented possi- 
bility of carrying out systematic investigations of the specific 
properties of each of these envirormients. An example of what 
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is now possible is demonstrated by the work of Povich et al. 
(2007). These authors used GLIMPSE and MSX (Midcourse 
Space Experiment, Price et al., 2001) data to perform a mul- 
tiwavelength study of the M17 complex, and found that the 
SED from the H II region peaks at shorter wavelengths and 
has a qualitatively different shape than the SED from the PDR. 
They also showed that different gas and dust components are 
mostly distributed in mutually exclusive regions: hot plasma 
occupies the very inner part of the H II region; warm ionized 
gas defines the cavity walls; warm dust permeats the ionized 
gas; and Poly cyclic Aromatic Hydrocarbons (PAHs) surround 
the whole star forming complex, tracing the PDR. 

In this paper, we make use of the newly available data to 
investigate the IR properties of the interior of evolved H II re- 
gions. We also address the question - from the observational 
point of view - whether dust grains survive the extreme con- 
ditions (i.e. intense UV photon flux and radiation pressure, 
stellar winds and high temperatures) which characterize this 
type of environments. Remarkably, evolved H II regions, due 
to their large physical size, are ideal candidates for carrying 
out such analysis, as they allow us to take fuU advantage of 
the high spatial resolution (e.g. from a few to tens of arcsecs) 
of these new data sets. 

The paper is organized as follows. In Section 2 and 3 we 
illustrate the data base and the source selection criteria. In 
Section 4, we generate source SEDs and describe the mod- 
elling procedure. In Section 5, we discuss our expectations, 
based on observations and theoretical predictions, of finding 
dust associated with a warm (lO'* K) gas. In Section 6 we 
show the relative spatial distribution of different populations 
of grains inferred by our multiwavelength data set. In Sec- 
tion 7 we compute the IR excess for each source and discuss 
its implications. In Section 8 we create color-magnitude and 
color-color plots, and show how these can be used for the 
identification of unresolved H II regions with respect to other 
classes of sources. We provide final remarks and conclusions 
in Section 9. 

2. THE DATA 

The bulk of the data is from the Herschel Hi-GAL survey. 
The Herschel infrared Galactic Plane Survey (Hi-GAL; Moli- 
nari et al. 2010) consists in the observations of the first and 
fourth Galactic quadrants using the PACS (Poglitsch et al. 
2010) and SPIRE (Griffin et al. 2010) instruments. The ob- 
servations were carried out in parallel mode at 70 fxm and 160 
fxm (PACS), and at 250 /xm, 350 /xm and 500 /xm (SPIRE), 
with an angular resolution from 6" (70 /im) to 35" (500 /im). 
During the Herschel Science Demonstration Phase (SDP), 
two Hi-GAL fields were observed. The fields are centered, 
for ^ = 0°, at / = 30° and / = 59°, and extend 2 degrees in 
both latitude and longitude. The data were processed with the 
ROMAGAL pipeline (Traficante et al., 2011), which allows 
for an accurate reconstruction of both compact and extended 
emission. Since the SDP data were taken early in the mission, 
corrective factors specific to extended emission were applied. 
In particular, PACS 70 /xm and 160 /xm data were divided by 
1 .05 and 1 .29, while SPIRE 250 fim, 350 /xm and 500 ^m data 
were multiplied, respectively, by 1.02, 1.05 and 0.94. For the 
calibration uncertainties, we assumed 20% for both PACS and 
15% for SPIRE (Bemard et al. 2010). 

We complemented the far-IR and sub-mm Herschel obser- 
vations, with MIPS 24 ^m (6") and IRAC 8 /im (2") data ob- 
tained by the MIPS Galactic Plane Survey (MIPSGAL) and 
by the Galactic Legacy Infrared Mid-Plane Survey Extraor- 



dinaire (GLIMPSE) survey. We adopted a 10% calibration 
error at both wavelengths (S. Carey, private communication). 
To overcome the occasional problem (e.g. the core of G39.8- 
0.3) posed by hard-saturated pixels in the MIPS 24 /xm data, 
we replaced the MIPS 24 /xm data with the combined MIPS 
24 /xm and MSX 21 /xm images produced by the MIPSGAL 
team and available to the members of the consortium. These 
images were generated taking into account appropriate color 
corrections and wavelength scaling (S. Carey, private commu- 
nication). 

We also made use of archival single-dish 6 cm data (see 
Section 3), which provided the basis for our source selection, 
and allowed us to constrain the evolutionary phase of the H II 
regions in the sample. In addition, we used the Multi-Array 
Galactic Plane Imaging Survey (MAGPIS) high-resolution 20 
cm data for a more accurate determination of the boundary of 
the H II regions and to perform our photometric measurements 
(Section 4). 

3. SOURCE SELECTION AND DEFINITION OF THE 
SAMPLE 

Since the launch of Spitzer in August 2003 (Werner et al. 
2004) and then of Herschel in May 2009 (Pilbratt et al. 2009), 
a steady flow of papers addressed multi-wavelength high- 
resolution investigations of Galactic H 11 regions. Apart for a 
few exceptions (e.g. the aforementioned Povich et al. (2007) 
work on M17, or the investigation of Flagey et al. (2011) of 
M16), all these studies (e.g. Watson et al., 2008; Deharveng 
et al. 2010; Zavagno et al. 2010; etc.) targeted bubble H n 
regions only. These sources are typically characterized by a 
common IR "structure": 8 /tm emission, dominated by PAH 
features, encloses the bubble; the bubble itself is filled with 
24 /xm emission, and this spatially correlates with free-free 
emission associated with ionized gas. In this work, one of our 
goals is to investigate if it is possible to extend and generalize 
to other classes of H II regions the results obtained from the 
analysis of bubbles. For this purpose, we focused our atten- 
tion on H II regions spanning a variety of morphological types. 
In order to assure uniformity of the sample, we required, for 
the selected sources, the availability of ancillary data allowing 
us to assess their evolutionary stage. 

The starting point of our selection was the catalog of Pal- 
adini et al. (2003). While for 29° < / < 31°, \b\ < 1°, the 
catalog lists 30 H II regions, for 58° < / < 60°, \b\ < 1°, we 
found only two sources. Therefore we concentrated our anal- 
ysis on the Hi-GAL SDP field centered at / = 30°. Out of 
the 30 cataloged sources falling in the 1 = 30° field, we made 
a further selection, by favoring those sources for which co- 
ordinates are known at least at the arcmin precision, and an 
estimate of the distance is available. We note that the distance 
information is central to our study, since it not only allows us 
to compute the luminosity of the sources, but also makes it 
possible to generate the linear longitude and latitude profiles 
discussed in Section 6. The final source selection consists of 
16 objects, including well studied H II regions such as N49 
(G288-0.2, Churchwell et al. 2006; Everett & Churchwell 
2010; Draine 2011) and W43 (G30.8-0.3, Bally et al. 2010). 
Table 1 gives details relative to each source, namely: source 
name (Column 1); Galactic coordinates (Column 2 and 3); 
6 cm flux (Column 4); Galactocentric distance (R, Column 
5); heliocentric distance (D, Column 6); far distance solu- 
tion in case of distance ambiguity (Dy, Column 7); electron 
temperature (T,,, Column 8); emission measure derived from 
single-dish (EMgcm. Column 9) and MAGPIS measurements 
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(EM20cm, Column 10); emission measure obtained assuming 
shell-like geometry of the source (EM,/,^;/, Column 11); elec- 
tron density computed from single-dish (n,, 6<;m> Column 12) 
and MAGPIS measurements (ne 20cm. Column 13); electron 
density computed assuming shell-like geometry of the source 
i^e.shdh Column 14). Table 2 provides additional information 
on angular sizes (9) and linear diameters (d). 

6 cm fluxes and angular diameters are from Altenhoff et 
al. (1979), Downes et al. (1980) and Kuchar & Clark 
(1997). Galactocentric and solar distances are, for aU the 
sources with the exception of three cases (G29. 0-0.6, G30.2- 

0. 1, G30.5H-0.4), from Anderson & Bania (2009, hereafter 
AB09), based on Lockman (1989) RRL observations. For 
these sources, AB09 solved the kinematic distance ambiguity, 
using HI emission/absorption and self-absorption data. We 
adopted their recommended values. In the case of G29.0-0.6, 
AB09 suggested that the source is situated at the far distance 
(11.5 kpc), although they obtained contradictory results from 
the HI emission/absorption data (near solution) with respect 
to the self-absorption data (far solution). In their paper, the 
authors stated that the HI emission/absorption measurements 
are usually more reliable to discriminate between near and far 
solution. However, for this specific source, the self-absorption 
data are of better quality, so they assigned it to the far distance. 
We visually inspected the '^CO data from the Galactic Ring 
Survey (GRS, Simon et al. 2001) and noted that G29.0-0.6 ap- 
pears to form a coherent structure with a companion source, 
G29. 1-0.7, for which the likely source of ionization (the Bl 
II type star S65-4) has a known photometric distance of 3.6 
kpc (Forbes 1989). For this reason, we re-located G29. 0-0.6 
to the near distance (3.4 kpc). For two sources (G30.2-0.1 and 
G30.5H-0.4) AB09 did not solve the distance ambiguity. For 
these, we computed the Galactocentric and solar distances us- 
ing the McClure-Griffiths & Dickey (2007) rotation model, 
taking Rq = 8.5 kpc and Oq = 220 km/sec, combined with the 
Lockman (1989) RRL measurements. The near and far solar 
distances for these two sources are quoted in Table 1 and both 
are used for the derivations of d, T^., EM and n^. For four 
sources the value of T^ reported in the table is from Downes 
et al. (1980 - G30.2-0.1, G30.6-0.1, G30.7-0.3) and Quireza 
et al. (2007 - G30.8-0.0). For the other sources, a crude esti- 
mate of Tg is obtained by applying the relation (Deharveng et 
al. 2000): 

r, = (372 ±38)/? + (4260 ±350) (1) 

For all the sources EM is computed from Mezger & Hen- 
derson (1967): 

Ti, = 8.235 X 10"2flr-'"^i^-2 '£M (2) 

In equation (2), T/, (the observed brightness temperature) 
and Te are in K, f is in GHz, and a is a constant set equal to 

1. EM, combined with the linear diameter, is used to derive 
Tie from: 

ne = ^^ (3) 

In the expression above, we assumed that n^. is uniform 
across the source and that this is spherically symmetric. How- 
ever, most of the sources in our sample do not have a spher- 
ical morphology. In Section 6 we will return on this point 
and dicuss its implications. Because EM has a weak depen- 
dence on Tg (EM oc T^^^), we do not expect the application 



of the empirical relation in equation (1) to bias significantly 
the derivation of this quantity and, in turn, of the electron den- 
sity. We note that all the emission measures quoted in Table 1 
are of the order of 10^ cm~^, suggesting that our sample is 
homogeneous and consists only of evolved H II regions. 

Finally, we computed, from the single-dish measurements, 
the Ly(Q:) luminosity of each source: 



LiLya) = 3.20x10 




(Garay et al. 1993) and the number of ionizing photons 
needed to excite the H II region: 

/ J' \ -0-45 

NLyc = 6.3 X lO^^photons s"M ^ j • (5) 

VGHz; 1027 erg s-iHz-1 

(Condon 1992). In the relations above, = Si^AttD^, and 
D and Tg are, respectively, the 6 cm flux, heliocentric dis- 
tance and electron temperatures which appear in Column 4, 
6/7 and 8 of Table 1. L(Lya) and Nl^c are reported in Ta- 
ble 6. The sources in our sample have ionizing luminosities 
in the range ~ 

1q48 _ iq50 photons s-\ with a mean value 
of 10'*^ photons s"'. According to Smith, Biermann & 
Mezger (1978), H II regions with an ionizing luminosity at 
least 4 times greater than the Orion Nebula (i.e. 4x 10"^^ pho- 
tons s~' at Rq = 10 kpc) are defined as giant. When prop- 
erly re-scaled by a factor 0.85 for Rq = 8.5 kpc (McKee & 
Williams 1997), the thereshold for giant H II regions becomes 
equal to 3.4x 10^^ photons s~', indicating that our sample con- 
sists of two giant H II regions (G29. 1+0.4 and G30.8-0.0) and 
14 sub-giant, many of which (at least 7) are close to the giant 
regime. 

4. SEDS 
4.1. Evaluating the SEDs 

The single-dish archival data used to establish the evolu- 
tionary stage of our sources have too coarse a spatial resolu- 
tion for allowing the accurate definition of the boundaries of 
the ionized-gas dominated region within each source. To this 
end, we made use instead the MAGPIS 20 cm data, which are 
characterized by an angular resolution of 6", i.e. compara- 
ble to the PACS 70 /im resolution. For each H II region, we 
drew MAGPIS contours at 2>-a above the background level 
and used these to evaluate the extent of the actual H II region. 

The IR wavelengths we are considering, as well as being 
contributed to by dust emission, are contaminated, especially 
in the PACS and SPIRE bands, by additional free-free (e.g. 
thermal bremstrahlung) emission. To evaluate this contribu- 
tion, we need to take into account the dependence on both 
frequency and electron temperature of the free-free spectral 
index, (S ^ A^^O- This is because the canonical approx- 
imation j3ff = 0.1 holds true only up to 10 GHz (3cm) and 
electron temperatures of the order of 8000 K. At higher fre- 
quencies (shorter wavelengths) and lower electron tempera- 
tures, the spectral index can differ from this value by a few 
to tens of a percent. Therefore, to each source we applied the 
relation (Bennett et al. 1992): 
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'^"''^ 11.68+ 1.51n(r^/8000/:)+ln(A/m) ' 

We used the values of Tg reported in Table 1 , together with 
the assumption that the abundance of ionized He (He++) is 
negligible with respect to the abundance of ionized H (H+). 
Table 3 provides, for each source, ^// in the range 24 /xm 
- 500 fira. We did not compute ^// at 8 fira, since in the 
mid-IR free-free transitions become sub-dominant with re- 
spect to free-bound transitions and these, in turn, are negli- 
gible when compared to dust emission (Beckert, Duschl & 
Mezger, 2000). From Table 3, we see that, in the range 24 iim 
to 500 /um, Pff varies between a minimum of 1.7 and a maxi- 
mum of 0.2, thus the subtraction of the free-free contribution 
becomes relatively important (of the order of a few percent) 
only at the longer SPIRE wavelengths. For each source, the 
free-free corrected flux at a given wavelength A, Sj, was ob- 
tained by subtracting from the IR flux, S^, the MAGPIS 20 
cm flux rescaled to A using: 

~ ^ ^ c /20xl0n''^^-^ 

Sx = Sj- 520cm* I -j^ I (7) 

A=20cm ^ ^ 

with A in jim. Noteworthily, the single-dish fluxes quoted 
in Table 1, properly rescaled to 20 cm using a canonical free- 
free spectral index 0.1, are within 20% of the MAGPIS 

measurements. 

After free-free subtraction, using the apertures obtained 
from the 20 cm MAGPIS data, we estimated the flux and built 
the SEDs by combining the MIPS 24 /im (and when, neces- 
sary, the MIPS/MSX data) with the PACS and SPIRE data. 
We also tried to include the 1.1 mm data, at 33" resolution, 
from the Bolocam Galactic Plane Survey (BGPS, Rosolowsky 
et al. 2010). However, the SED modelling revealed a sig- 
nifcant flux loss at this wavelength (Figure 1 and Figure 2). 
The origin of this effect can be attributed to the BGPS data 
processing, which largely (at more than 90% level) removes 
emission on scales of or greater than 5.9'. This operation is 
necessary in order to separate the astronomical signal from 
the atmospheric fluctuations in the data streams. For struc- 
tures larger than 3.8', the attenuation is of the order of 50% 
(Aguirre et al. 201 1). The 16 H II regions in our sample have 
an average angular size of 3.2', hence still in the range in 
which the extended emission should be preserved for the most 
part. Therefore, our result appears to indicate that the amount 
of filtered flux is more than expected, at least on the scales we 
are considering. We could compensate for this effect by arti- 
ficially increasing the systematic error. Instead, we opted for 
a more conservative approach and decided to exclude the 1.1 
mm data from the fit. We note that the loss of signal on large 
angular scale would have been comparable or even greater if, 
rather than the BGPS data set, we had used the 19" resolu- 
tion ATLASGAL (Schuller et al. 2009) data at 870 /^m. In 
this case, a significant fraction of flux is removed from the 
data already on scales as small as 2.5', that is smaller than the 
average size of the H II regions in our sample. 

For 24 /zm < A < 500 /im, the flux at each wavelength 
and for each H II region was determined within the defined 
apertures after point source subtraction - performed with he 
Starfinder algorithm (Diolati et al., 2000) - and convolution 
to the SPIRE 500 angular resolution (35"). To evaluate 
the background level, we took the median within a control re- 
gion, selected to be close to the source under consideration. 



yet paying attention to avoid potentially contaminating bright 
objects. Uncertanties on measured fluxes ((Ts), prior to free- 
free subttaction, were obtained by adding in quadrature the 
cahbration error and the background standard deviation, ac- 
cording to: 

/ N \ 1/2 

os=\^^sSn'oI^^ (8) 

where 5, is the flux in pixel i, N is the number of pixels 
within the aperture, and back is the median background in 
the control region. To this error, due to the subtraction of the 
free-free contribution, we added in quadrature an extra 30% 
uncertainty associated with the MAGPIS data (e.g. Povich et 
al. (2007)). We also removed the zodiacal light contribution 
from the 24 /im data. This operation was carried out using the 
entry 'ZODY_EST' in the FITS header of each 24 /im frame, 
where the amount of zodiacal light per frame was estimated 
by the Spitzer Science Center from the predictions on Kelsall 
et al. (1998). 

The photometric measurements performed after back- 
ground subtraction retrieved, at every wavelength and for each 
source, a positive flux (see Table 4). The average S/N (Signal- 
to-Noise) across aU the bands and for the entire sample is 
S/Nx^,i = 3.6±2.1. If we break this number per wavelength, 
we have: S/Nsf,„ = 4± 2.3, S/N2A^,m = 6.7± 1.9, S /Njo^.^ = 
3.3± 0.9, S/Nmjun = 2A± 1, S/N25q^„, = 2.9± 1.6, ^/A^jso^m 
= 2.6± 1.4 and S/Nsm^im = 2.7± 1.5. We note that, as a gen- 
eral trend, the significance of the detection tends to decrese 
with increasing wavelength. This is hkely due to the fact that 
at A ~ 160 /im, the emission from the H II region starts com- 
peting with the emission of the Interstallar Medium (ISM). At 
longer wavelengths (i.e. in SPIRE bands), this contribution 
becomes sub-dominant with respect to the local background. 

It is worth mentioning that the detected IR emission that we 
ascribed entirely to H II regions could instead originate from 
warm/cold material located in the foreground/background 
along the line of sight. Indeed this bias at least partly affects 
our flux extraction procedure. In fact, in the presence of a 
structured, highly- varying background such as the one char- 
acterizing the Galactic Plane, background- subtraction does 
not assure the complete removal of emission unrelated to the 
source, due to the intrinsic difficulty of estimating the back- 
ground level. However, what provides confidence in our re- 
sults is the fact that, for the majority of the sources, the MAG- 
PIS 20 cm emission, which traces the ionized gas component 
of our H II regions, does not (or not completely) overlap with 
the peak of the IR emission in the other bands (likely asso- 
ciated with the PDR), suggesting that the H II region and its 
dust content is largely (or partly) exposed to the observer. 

4.2. SED Fitting 

Figure 1 and 2 show the SEDs derived from the photometric 
measurements described in the previous section. Since the pi- 
oneering works of Chini et al. (1986a, 1986b, 1986c, 1987), it 
is known that, in order to explain observations of H II regions 
above and below ^100 /im in the framework of traditional 
dust models, that is with a constant spectral emissivity index 
I3dust< it is necessary to invoke the existence of a 2-temperature 
component dust distribution: a warm, low density population 
of dust grains situated in the proximity of the ionizing source, 
and a colder dust population far from the central star (or stars). 
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One has to keep in mind that the observations that Chini et 
al. used to constrain their modelHng were taken in the early 
Eighties and, as such, were characterized by relatively low 
angular resolution, with the consequence that the actu£il H n 
region was not resolved with respect to the associated PDR. 

Both laboratory experiments and observations conducted in 
recent years suggest that traditional dust models assuming a 
constant spectral emissivity index might need to be revised 
(see discussion below). For this reason, we tried to fit our 
SEDs in three possible ways: (1) with an isothermal modified 
blackbody with constant spectral emissivity index (Pdusi = 2); 

(2) with an isothermal modified blackbody with variable Pdus, ; 

(3) with a 2-temperature component model with (3^ust =2- We 
did not attempt the fit with a 2-temperature component model 
with variable ^dust since the measurements would not provide 
enough constrains (6 data points) with respect to the number 
of free parameters (6 parameters). We also did not include 
in the fit the IRAC 8 /im band, which is typically contributed 
to by emission associated with stochastically heated R\Hs. 
To properly take into account PAH emission, we would need, 
for each H II region in our sample, an accurate knowledge of 
the radiation field generated by the central source, and this 
information should be coupled to a sophisticated dust model 
(e.g. DustEM, Compiegne et al. 2011) able to reproduce, in 
a coherent way, both the aromatic features and the continuum 
thermal emission due to larger grains. Since detailed informa- 
tion on the radiation field is not available for most of our H II 
regions, we considered only the data at wavelengths longer 
than 24 fim. The three functional forms corresponding to our 
three adopted models are: 

sZ=Ai,i(^^^ Bx(Tdus,A) (9) 
for an isothermal model with constant emissivity index; 
/ A \ 

5A.2=Ai,2f^j BxiTdus,,!) (10) 

for an isothermal model with variable Pdust, and: 

S^3=Ai,3(^^^ Ba(7;,3)+^2,3 BxiT„,3) (11) 

for a 2-temperature component model, where T^.^ and Tj^ 3 
are, respectively, the temperatures of the cold and warm com- 
ponents, and Ao is set to 100 /im. 

The choice of iS^ust = 2 is motivated by the fact that ob- 
servations of the silicate absorption feature near 10 iim seem 
to indicate that BGs are amorphous materials (Kemper et al. 
2004), and these are expected to behave in a similar fashion 
to crystalline dielectric materials, for which sub-mm absorp- 
tion has a temperature-independent quadratic dependence on 
frequency (Boudet et al. 2005). 

To estimate the parameters Ai_i, Tj^st,!, M,2, Tdust,2, M,3, 
Tc,3 and T^ s for each model (equation (9), (10) and (11)), we 
used a Monte Carlo Markov Chain (MCMC) method (Lewis 
& Bridle, 2002). This technique consists in looking for the 
likelihood maximum by sampling the parameters space with a 
Metropohs-Hastings algorithm. The a priori probability den- 
sities of the parameters were set to be as wide as possible to 
avoid introducing artificial biases in the fit. Color corrections 
were computed iteratively during the fitting procedure using 
the MIPS, PACS and SPIRE transmission filters. In addition. 



systematic (i.e. calibration) uncertainties were taken into ac- 
count and included in the parameters error bars by means of a 
dedicated Monte Carlo run. A straightforward goodness- 
of-fit method was applied in parallel with the MCMC method 
to test the quality of the fits. In general, the best-fit values re- 
covered by the two methods (MCMC and goodness-of-fit) 
were found in agreement within a few %. 

The MCMC results of the fitting procedure for the isother- 
mal model with varying spectral emissivity index and the 2- 
temperature component model are given in Table 5. We do not 
show the results of the isothermal model with fixed f3dust as 
this model fails to reproduce simultaneously the 24 /vm - 500 
/im observations. For the remaining two models, we made 
trial fits by including/excluding the 1.1 mm BGPS measure- 
ments. In both cases, the best-fit parameters (i.e. the best x^) 
were obtained by excluding the 1.1 mm Bolocam extracted 
flux, as discussed in the previous section. A direct comparison 
between the isothermal and 2-temperature component models 
in terms of could not be perfomed, since the degrees of 
freedom are not constant. However, a comparison in terms 
of the Cumulative Density Functions (CDF)" for the x^ dis- 
tribution of each model model revealed, for all the sources 
but three, that the 2-temperature component model is signifi- 
cantly more in agreement with the observations (see Column 
7 and 13 in Table 5). For G30. 1-0.2, the 1-temperature model 
seems to work surprisingly better than the 2-temperature one, 
although the best-fit fi^ust (0.32) lies outside the normally ac- 
cepted range of values. In the case of G30. 3-0.0 and G30.6- 
0.1 none of the models appears to be a good representation 
of the measured SEDs. For the 13 sources for which the 2- 
temperature dust component model is in agreement with the 
data, we obtain that the cold component peaks between 100 
/xm and 160 /im and is characterized by a temperature in the 
range ~ 20 - 30 K, while the warm component has a peak 
around 20 /xm - 70 /xm with a temperature from ~ 50 to ~ 
90 K. The average values for the cold an d wa rm components 
are, respectively, 7;,3 = 24.53±5.1 K, and r„,3 = 68.43±20.52 
K. We note that T^j does not differ significantly from temper- 
atures measured in PDRs associated to WBB (e.g. Rodon et 
al. 2010, Anderson et al. 2010). In principle, both cold and 
warm components could be attributed to BGs. This is indeed 
possible only if the warm and cold dust emissions are not spa- 
tially correlated, since a scenario in which two populations of 
large grains, subject to the same radiative field, thermalize at 
very different temperature is rather unreaUstic. An alternative 
solution, for spatially coincident components, would be to as- 
sign the cold component to a population of BGs, and the warm 
component to a population of VSGs. This option would be in 
agreement with the increase of VSG abundance in ionized en- 
vironments reported (see Section 5) by Paradis et al. (2011) 
and Flagey et al. (2011), and which can be explained through 
partial sputtering of the BGs in the most inner parts of the H II 
regions. Without a more refined dust modeling which can take 
into account the radiation field properties of the sources, we 
cannot unambiguously conclude on this point. 

For the 2-temperature component model, we used a con- 
stant Pdusi- However, an increasing amount of evidence ob- 
tained both from laboratory experiments (Agladze et al. 1996, 
Mennella et al. 1998, Boudet et al. 2005) and observations 
targeting the Cosmic Microwave Background (Dupac et al. 

" The CDF (or 1-p-value) gives the probability that a random variable X 
with a given probabiUty distribution will be found at a value less or equal than 

X. 
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Fig. 1. — SEDs for all the H II regions in our sample. Triangles denote photometric measurements at 24 fim (MIPS), 70 fim and 160 fim (PACS), 250 fim, 
350 fj,m and 500 fim (SPIRE). Red squares indicate the IRAC 8 fj,m and Bolocam 1.1 mm data points. Also shown are photometric uncertainties. The solid line 
corresponds to the best-fit isothermal model with varying emissivity spectral index. 
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2003, Desert et al. 2008, Paradis et al. 2010; Veneziani et al. 
2010) appears to support a temperature dependence (as well 
as a wavelength dependence, e.g. Paradis et al. 201 lb) of the 
spectral emissivity index. Theoretically, the existence of an 
inverse relation between dust temperature and pdust could be 
interpreted in light of the model recently proposed by Meny 
et al. (2007). This model, which was built over the Two Level 
System (TLS) theory first formulated by Phillips (1972) and 
Anderson et al. (1972), explains the observed T^^st - Pdusi 
anti-correlation as due to both the disordered charge distri- 
bution characterizing amorphous materials and TLS (or tun- 
neling) effects which occur on quantum scale in the grains 
structure. 

To analyze the possibility of a variation of (3dust with tem- 
perature, we considered the cold component only, which we 
fitted with a modified blackbody. For the warm component, 
we do not have enough measurements to constrain the fit with 
three free parameters i(3dust, ^2,3, Twj)- Since for the ma- 
jority of the sources (10 out of 16) the 70 /zm flux seems to 
contribute to the cold component (Figure 2), we included all 
the measurements in the range 70 ^m < A < 500 fim. 

The major complication in investigating a potential inverse 
relation between the spectral emissivity index and dust tem- 
perature is that these two quantities are intrinsically degener- 
ate in parameter space. In fact, their 2-D posterior probabil- 
ity has an elongated, slant shape, and is characterzed by the 
same functional form of the anti-correlation that we are try- 
ing to identify in the data (Shetty et al. 2009). In order to 
overcome this issue, we adopted the following procedure. We 
first fitted the measured fluxes using the MCMC algorithm, 
which allowed us to recover the joint posterior distribution of 
the parameters, thus keeping under control their intrinsic de- 
generacy. We then applied a technique successfully tested on 
BOOMERanG data (Veneziani et al. 2010), consisting in es- 
timating the Tdust - l^dust relation through a Monte Carlo of 
the caUbration errors. That is, for each iteration 7 (100 in to- 
tal): (1) we fitted equation (10) to each source and derived the 
best-fit values of the parameters; (2) we then estimated the 
relation: 

^.«.=Ax(||)" (12) 

by combining the results obtained for all the sources; (3) at the 
end of the 100'* iteration, we had a set of 100 pairs of Aj and 
aj values, one for each realization of the calibration error. By 
marginalizing over the calibration errors'^, we obtained the fi- 
nal best guess for A and a: A = 3.0 ± 0.2, a = -1.5 ± 0.2. At 
first sight, this result might seem to confirm the existence of 
an inverse relation between spectral emissivity index and tem- 
perature. However, several authors (Masi et al. 1995, Shetty 
et al. 2009) argued that the observed anti-correlation might be 
spurious and due to Une-of-sight temperature mixing effects. 
To test this hypothesis, we performed a simple simulation. We 
generated a synthetic random distribution of isothermal (30K) 
H II regions and ISM, for which we adopted a dust temper- 
ature of 17 K (Boulanger & Perault, 1988b). For both the 
H n regions and the ISM, we took pdust =2. With the same 
strategy described above, consisting in applying the MCMC 
technique for both fitting the sources SEDs and estimating the 
best-fit values A and a, we obtained: A = 3.1 ±0.1 , q: = -1.1 

This operation is performed using the GetDist software of the public 
CosmoMC package. 



± 0. 1. These values are very similar to the ones derived from 
the 16 evolved H II regions. If, instead of /Sdusi = 2, we adopt a 
shghtly lower value, such as Pdust = 1-7, then the best-fit anti- 
correlation for the simulated sample almost perfectly ovelarps 
with the inverse relation estimated from the real sample (see 
Figure 3). We interpret this results as a strong indication that 
the observed anti-correlation is, at least in this case, Ukely due 
to a line-of-sight temperature mixing effect. 

5. CAN DUST SURVIVE IN H II REGIONS ? 

The presence of dust in the warm ionized gas has been ex- 
tensively debated in the literature, due to its relevance for our 
understanding of the processes of destruction and formation 
of dust grains in the ISM. Given the extreme conditions char- 
acterizing the interior of H II regions, theoretical models pre- 
dict the depletion or partial destruction of the grains in these 
environments. Three main mechanisms could be responsible 
for this phenomenon (Inoue 2002): (1) radiation pressure; (2) 
stellar wind by the central source; (3) dust sublimation. The 
effect of radiation pressure was investigated by Gail & Sedl- 
mayr (1979b), and more recently by Draine (2011). Accord- 
ing to Gail & Sedlmayr (1979b), radiation pressure is capable 
of generating a central zone of dust low-density (or cavity) of 
size of the order of 20% the ionization radius. Gas cavities 
are instead known to be produced by stellar winds (Comeron 
1997). However, given the coupling between dust and gas, a 
corresponding dust depletion zone is also expected. Sublima- 
tion appears to be the least effective process. In fact, if on one 
side, subhmation of dust grains could indeed occur as a con- 
sequence of the intense radiation field characterizing the most 
inner parts of an H II region, yet the estimated radius of the 
cavity induced by this mechanism is only of the order of 10^ 
pc (Mookerjea & Ghosh 1999). 

From the observational point of view, one of the main find- 
ings of Povich et al. (2007) is the significant decrease of PAHs 
inside M17, which is thought to be due to destruction of the 
aromatic molecules due to the extreme ultraviolet (EUV) flux. 
Indication of PAH depletion associated with ionized gas is 
also reported by Paradis et al. (201 1) for the Large Magellanic 
Cloud. These authors also found that the relative abundance 
of VSGs with respect to BGs appears to increase when transi- 
tioning from the neutral to the ionized medium. Flagey et al. 
(201 1) arrived at a similar conclusion by analyzing the bright 
24 /im emission filling up the interior of the M16 H II re- 
gion, which they argued was caused by sputtering of BGs into 
VSGs. Finally, while the existence of BGs in neutral PDRs is 
supported by several studies (e.g. Povich et al. 2007; Com- 
piegne et al. 2007), their survival in the interior of H II region 
is rather controversial and so far it is corroborated only by 
statistical analysis, such as cross-correlation of large data set 
tracing both IR and free-free emission (Lagache et al. 2000; 
Paladini et al. 2007, Planck Collaboration, 2010). 

6. IR EMISSION DISTRIBUTION INSIDE AND 
AROUND H II REGIONS 

In this section, we address the question of if stnd what kind 
of dust is present in H II regions, and investigate the interplay 
between the neutral PDR and the inner H II region. Follow- 
ing the same procedure outlined in Section 4, we used the 
MAGPIS 20 cm contours to trace the ionized gas content of 
the inner H II region, and IRAC 8 /xm, MIPS 24 /xm, PACS 
70 /im, and SPIRE 250 data to trace PAHs, VSGs and 
BGs potentially associated with both the PDR and the H II 
region. With this radio/IR data set, we created two types of 
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Fig. 2. — SEDs for all the H II regions in our sample. Triangles denote photometiic measurements at 24 fim (MIPS), 70 /im and 160 fim (PACS), 250 fim, 350 
fim and 500 /im (SPIRE). Also shown are the IRAC 8 fj,m and Bolocam 1.1 mm data points (red squares) and photometric uncertainties. The solid line indicates 
the best-fit 2-temperature dust component model with fijust =2. Dashed lines illustrate the warm and cold components. 
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Fig. 3. — Simulated sample of a random distribution of H II regions and 
ISM (grey points and error bars) overlaid on the observed H II region sample 
(red points and error bars). The grey solid line is the best-fit function for the 
synthetic sample, while the red dashed line illustrates the best-fit model for 
the 16 H II regions considered in this work. 



3-color images (Appendix A). The first set of images (left 
panels in Figure Al) is obtained by combining MAGPIS 20 
cm (red), MIPS 24 fim (green) and IRAC 8 fim (blue) data, 
while the second set (right panels in Figure Al) is generated 
from the composite MAPGPIS 20 cm (red), SPIRE 250 fim 
(green), PACS 70 /im (blue) data. A visual inspection of these 
figures reveals a wide range of morphologies of the sample 
H II regions, namely: bubbles (e.g. G28.8-0.2, G29.0-0.6, 
G31.1H-0.3); elongated structures (e.g. G29.1H-0.4, G30.5- 
0.3); diffuse nebulae (e.g. G30. 1-0.2, G30.6-0.1); and more 
complex shapes (e.g. G30.8-0.0, G31.0H-0.5). 

Despite the complexity posed by dealing with different ge- 
ometries, we identified 3 general trends in the color images: 
(1) the radio and 24 /im emission ovelap in several cases; (2) 
the 8 /im and 70 fim emission are spatially correlated; (3) 
the 250 fj,m emission distribution is frequently confused with 
the background level. In regions where the ISM emission is 
more attenuated, the peak of the 250 /im emission is often dis- 
placed with respect to the 20 cm peak, although one can find 
instances of partial overlap. 

To investigate even further the relation between dust and 
ionized gas, we generated latitude and longitude emission 
profiles. For this purpose, we used MAGPIS 20 cm, IRAC 
8 fim, MIPS 24 fim and SPIRE 250 fim data. The profiles 
were obtained following these guidelines: (1) for each source 
and at each wavelength, we created a postage-stamp image, 
either 10' x 10', 15' x 15' or 20' x20', depending on the angu- 
lar extent of the source; (2) we convolved the postage-stamp 
images at the same angular resolution (18"); (3) we regridded 
the images using a 6" pixel size; (4) we generated the actual 
profiles by computing the median along each column (lati- 
tude cut) or row (longitude cut) of the convolved and regrid- 
ded postage-stamp images; (5) at each wavelength we nor- 
malized each latitude/longitude profile to its maximum; (6) 
finally, for every source, we overlaid the profiles at different 
wavelengths. The resulting profiles (see Appendix B) high- 
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Fig. 4. — Top to bottom panel: distributions of the linear separation be- 
tween the peaks of emission of, respectively, SPIRE 250 fim, IRAC 8 ^m, 
MIPS 24 fim with respect to MAGPIS 20 cm. 



light the relative spatial distributions of the various compo- 
nents of dust emission already evidenced by the 3-color im- 
ages: the 8 fim and 250 fim emission peak almost exclusively 
outside the gas-dominated region, while the prominent emis- 
sion at 24 fim appears to be closely associated with the H II 
region. To quantify this effect, we used the latitude and lon- 
gitude profiles to measure the linear separation between the 
peaks of emission at different wavelengths. To this end, we 
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excluded all ambiguous cases, i.e. H II regions for which the 
identification of the source with respect to the background, ei- 
ther in the latitude or longitude profile, is not straightforward. 
This is the case of G29.1-0.0, G29.2-0.0 and G31. 1+0.3 (lati- 
tude profile), G30. 1-0.2 (both latitude and longitude profiles), 
G30.5+0.4 (longitude profile). For the remaining sources, we 
computed from the profiles the angular distance of the IRAC 
8 /im/SPlRE 250 /im/MlPS 24 /im peak of emission with re- 
spect tp the MAGPIS 20 cm peak. With an estimate of the an- 
gular separation for each pair of wavelengths, we took the av- 
erage of the values obtained from all the latitude and longitude 
profiles, and used the distance of the source to convert such an 
average into a linear distance (in pc). Figure 4 shows that the 
average linear separation between the peak of the 20 cm emis- 
sion with respect to the peak of the 24 emission is half the 
distance between the peak of the same radio continuum emis- 
sion and the 8 /im emission and this, in turn, is also half the 
separation between the MAGPIS peak and the 250 /im one. If 
we assume, as customary, that the 8 /im data trace PAH emis- 
sion and the 250 /im data trace BG emission, we can read this 
result as a strong indication that PAHs and BGs are not spa- 
tially correlated, with the PAHs statistically distributed along 
a ridge which appears to be closer to the ionizing source than 
the BGs ridge. This stratification is consistent with the sce- 
nario recently described in Draine (201 1), in which radiation- 
pressure-driven drift is very effective in moving grains with a 
size a > 0.01 /im (BGs) outward, while smaller grains (PAHs) 
also drift but more slowly. Remarkably, this is also in agree- 
ment with Krumholz & Matzer (2009) who claim that radia- 
tion pressure is usually unimportant for H II regions ionized 
by a small number of stars but becomes a dominant factor for 
the expansion dynamics of large H II regions. What remains 
as very puzzling is the source of the 24 jim emission. If this 
wavelength traced a population of VSGs co-eval with respect 
to the BGs and PAHs at the periphery of the H II regions, we 
would expect from the predictions of Draine (2011) to find 
the emission peak between the 8 /im and the 250 /im peak. 
Likewise if the 24 /im emission was due to a co-eval distribu- 
tion of warm BGs. We then tentatively attribute this emission 
to a new generation of either VSGs or BGs. We note that 
this conclusion is similar to what was advocated by Everett & 
Churchwell (2010) from the analysis of the behaviour of the 
24 /im emission in WBB: that is, within bubbles, dust grains 
are re-supplied by destruction of embedded, dense cloudlets 
of interstellar material that are overrun by the expansion of 
the bubble itself. 

The latitude and longitude profiles also evidence that some 
of the sample H II regions are very prominent in emission 
at all wavelengths, while others are barely visible above the 
background level. This might indicate that different evolu- 
tionary stages are present in our sample. Although all the 
H II regions we are considering are evolved, some might be 
more evolved than others and, in particular, close to the stage 
where the internal pressure equalizes the external one, like, 
for instance, G30. 1-0.2. Alternatively, as suggested by the 
recent MHD simulations by Arthur et al. (201 1), highly irreg- 
ular morphologies might be associated with highly ionized, 
density-bounded nebulae, powered by the hottest stars. Cer- 
tainly, a mixture of both scenarios is also a possibility. 

Both the 3-color images and the profiles clearly show that 
the vast majority of the H II regions we are analyzing are not 
characterized by spherical symmetry. Therefore, the angu- 
lar diameters (and corresponding linear diameters) obtained 
from gaussian fits of the single-dish data (Column 2 of Ta- 



ble 2) have to be considered only as a rough indication of 
the actual size of the sources. We then used the MAGPIS 
20 cm contours to measure the angular size of each source 
along its major ((?,„„/) and minor (6',„,„) axis (Colum 4 and 5 
of Table 2). The analysis of the latitude and longitude pro- 
files also reveals that the 20 cm emission of 6 H II regions 
(G29.0-0.6, G29.1+0.4, G30.3-0.2, G30.8-0.0, G31. 1+0.5, 
G3 1.5+0.3) is characterized by a double-peak behavior, sug- 
gesting that these objects might present a shell-like geometry. 
If this was the case, the electron densities derived by assum- 
ing that a much larger volume is occupied by the ionized gas 
would be highly underestimated. Therefore, we computed the 
electron densities from the MAGPIS measurements both in 
the hypothesis of uniformly filled sources (case A) and, by 
assuming a shell-like structure for the 6 H II regions showing 
a double-peak profile (case B). For case A, electron densities 
were obtained using the geometric mean, c/^'^, of the linear di- 
ameters di„ai and ii,„,„ inferred from the MAGPIS angular ma- 
jor and minor axis. For case B, we estimated the angular size 
of the shell from one of the peaks of the 20 cm double-peak 
profile (Column 8 of Table 2) and we converted this value into 
a shell thickness {dsheii, see Column 9 in Table 2) using the 
heliocentric distance of the source. To evaluate the electron 
density, we doubled the estimated dsheii, in order to account 
for the double crossing of the shell with respect to the line of 
sight. The electron densities obtained for case A and B are 
provided in Table 1 . 
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Fig . 5 . — Electron density-size relation derived from: ( 1 ) single-dish 6 cm 
data (green circles, solid line), (2) MAGPIS 20 cm data (blue circles, dashed- 
dotted line), (3) MAGPIS 20 cm data (red circles, dashed line) taking into 
account the potential shell-type morphology of 6 sample HIT regions. 

With the derived electron densities and linear sizes, we ex- 
plored the n^-d relation. For a uniform-density nebula, n^, is 
expected to scale with the linear diameter as d~^/^. Depar- 
tures from this behavior can be explained by either invok- 
ing a non-uniform, clumpy nature of H II regions (Garay & 
Lizano, 1999; Kim & Koo 2001), or the presence of dust par- 
ticles in the ionize gas, which would absorb part of the ion- 
izing radiation and would cause the H II region to be smaller. 
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Here, for completeness, we considered the values of and 
d estimated both from the single-dish and MAGPIS data, and 
from both case A and B. Figure 5 shows the best- fit model 
for, respectively, (1) the single-dish measurements assuming 
that the ionized gas uniformly fiUs the source volume (green 
dots); (2) the MAGPIS measurements with the same approx- 
imation as in (1) (blue dots); (3) the MAGPIS measurements 
conidering the potential shell-like geometry of 6 sample H II 
regions (red dots). Visibly, when the shell-like geometry is 
adopted, the inferred electron densities (linear diameters) are 
~ one hundred times higher (~ ten times smaller) than the 
values obtained using correspondingly larger volumes. In the 
absence of these large/small electron densities/linear diame- 
ters, the Tie-d distribution is flatter and can be well represented 
by a power-law with spectral index > -1. However, when the 
effective size of the source is correctly evaluated, the best- 
fit relation (n^ ^ d"'^^) steepens significantly, ruling out the 
possibility that a significant amount of dust is present in the 
sample H n regions. 

7. IR EXCESS 

The presence of dust inside H II regions can also be inves- 
tigated in terms of energy budget, through the so-caUed In- 
frared Excess (Garay et al. 1993): 



IRE- 



L{Lya) 



(13) 



where L/^ is the IR luminosity and L(Lya) is the Lyman 
photon luminosity derived in Section 3. Evaluation of the 
IRE provides an estimate of the fraction of photons emitted 
by the central source and directly absorbed by dust particles. 
We computed the IRE both for H II regions only, and for the 
combined H II regions and associated PDRs. For H II regions, 
the IR luminosity was obtained by integration of the best-fit 2- 
temperature component models (5a,2) derived in Section 4.2: 



'^IRJIII ■ 



--4ttD' 



S\2 dX 



(14) 



In the expression above, we set Am,„ = 1 /xm and \„ax = 
lO'^jim, and is, for each source, the solar distance pro- 
vided in Table 1 . In order to estimate the IR luminosities for 
the (H II region + PDR) ensembles, we perfomed photomet- 
ric measurements between 24 /xm and 500 jjbm by including 
in the aperture not only the region identified by the MAG- 
PIS 20 cm contours, but also the PDR. For this purpose, in- 
stead of using the MAGPIS 20 cm data to generate our aper- 
tures, we used the IRAC 8/im images. Following the same 
steps described in Section 4, we estimated and removed the 
free-free contribution in each band, and fitted the resulting 
SEDs with a 2-temperature component model with jijust = 
2. LjR Hii+pDR for each source was then evaluated by apply- 
ing a relation similar to equation (14). The derived Ljr hii, 
^iR,Hii+PDR, IREfi// and IREhh+pdr are shown in Table 6. For 
all the sources (with one exception, i.e. G30.2-0.1), both 
IRE//// and IREhh+pdr are slightly (but consistently) greater 
than 1. However, IREhh+pdr is systematically higher than 
IRE////, with IREhii = 1-08 ± 0.14 and IRE,iii+pdr = 1.17 ± 
0.04. As expected, this result indicates that the amount of 
radiation absorbed in the PDR is larger than the amount ab- 
sorbed in situ. It also shows that, if dust co-exists with the 
ionized gas in the interior of H II regions, its mass is not sig- 
nificant, in agreement with what we found from the analysis 



of the Ue-d relation described in the previous section. Finally, 
we note that the values of IRE obtained for the present sample 
of H II regions are substantially lower (a factor 10 on average) 
than the IR excess reported by Garay et al. (1993). This is 
likely a consequence of the different evolutionary stage of the 
sources in the two samples: the sources we are considering 
are more evolved and less dusty, while the compact H II re- 
gions analyzed by Garay et al. (1993) are still embedded in 
the natal cocoon. 

The photometric measurements for both the H II regions 
and the (H II region + PDR) systems allowed us, as well as 
to derive IR luminosities, to carry out a qualitative compari- 
son between the average behaviour of their SEDs in the wave- 
length range 24 /im < A < 500 /xm. The average SEDs, for the 
H II regions and (H 11 region -I- PDR) systems separately, were 
obtained by first normalizing each individual SED to its peak 
and then by taking the average of these normalized SEDs. The 
result is illustrated in Figure 6. The peak of the H II regions 
SED is shifted towards shorter wavelengths (^ 70 /im) with 
respect to the (H II regions -i- PDRs) SED peak (~ 160 /xm). 
This trend is in agreement with the finding of Povich et al. 
(2007) and can be ascribed to the higher temperature of the 
dust grain population in the H II region compared to the tem- 
perature of the grain population located in the PDR. 

8. COLOR-COLOR PLOTS 

We used the fluxes extracted in the previous section for the 
(H II region + PDR) complexes to generate color-magnitude 
and color-color plots. Goal of this exercise was to look for 
characteristic colors which can allow the identification of this 
class of objects with respect to other populations of Galactic 
sources. In particular, we considered sources in either an ear- 
lier or later evolutionary stage with respect to our sample H II 
regions. In the former category (e.g. younger sources) fall 
the approximately 100 sources listed in the Hi-GAL prelimi- 
nary catalog (Elia et al., 2010) and located in our 2° x2° SDP 
field. These sources have measured fluxes at all PACS and 
SPIRE wavelengths and at 24 /im. Most of them are thought 
to be Young Stellar Objects (YSOs) (EUa et al. 2010). For the 
latter category (e.g. older sources) we made use of the 39 Su- 
pernova Remnants (SNRs) cataloged by Pinheiro-Gongalves 
et al. (201 1, hereafter PGll), based on GLIMPSE and MIPS- 
GAL data. 

The color-magnitude and color-color diagrams are shown 
in Figure 7. Panels b), c) and d) highlight the general segre- 
gation of the (H II regions + PDRs) with respect to the YSO 
population. This effect can be attributed to the relatively high 
dust temperatures (~ 30 - 50 K) characterizing the (H II re- 
gion + PDR) systems, which causes the peak of their SED to 
occur at wavelengths shorter than ~ 160 fim (see Figure 6 and 
Section 7). Conversely, in the very early stages of star forma- 
tion (e.g. YSOs) typical dust temperatures do not go beyond 
~ 15 K, with the consequence that the SED peak is located 
at A > 160 fim. In practise, while the SED of an (H II region 
+ PDR) is rising fast between 24 /xm and 70 /xm and drop- 
ping rapidly between 160 /xm and 250 /xm (and between 250 
/xm and 500/im), that of a YSO is also rising in the 24 /xm 
< A < 70 /im but much more slowly, while for 160 /im < A < 
250 fim (and 250 /im < A < 500) it keeps steepening or, if it 
drops, it drops slowly. Since the Hi-GAL preUminary catalog 
was not band-merged with the GLIMPSE 8 /xm catalog, we 
were not able to compare the mid-IR color distribution of the 
(H II region + PDR) population with that of YSOs. 

Panel a) illustrates the behaviour of the (H II region -t- PDR) 
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systems compared to that of SNRs in the 8 /im/24 /im vs. 
24 jimllQ jiva color-color space. Both populations appear to 
span quite a wide range of 8 jimllA /zm ratios. However, 
while SNRs also strech along the x-axis direction for an or- 
der of magnitude, the (H II regions -I- PDRs) are confined to a 
narrow range of 24 /xm < A < 70 values. The relatively 
spread out distribution of SNRs likely reflects the variety of 
morphologies and emission mechanisms highhghted in PGl 1. 
It is worth to mention that a partial overlap of mid-IR colors 
for the composite (H II region + PDR) systems and SNRs was 
already reported by Arendt (1989). 

A search for Hi -GAL counterparts of the SNRs in the PGl 1 
catalog is currently on-going (Noriega-Crespo et al., 2012 in 
preparation). Preliminary results suggest that the detection 
rate in the PACS 70 ^im band should match the MIPSGAL 
one, while only few, very bright SNRs (e.g Gil) are detected 
at 160 /im, and almost none in the SPIRE bands. This ef- 
fect is currently interpreted as either due to an intrinsic lack 
of emission from these sources at far-lR wavelengths or to 
a competing background level (Noriega-Crespo, private com- 
munication). 




[/im] 

Fig. 6. — Average H 11 region (solid line) and PDR (dashed line) SED. The 
average is taken over the sample of 16 sources considered for this analysis. 



9. CONCLUSIONS 

The analysis of Spitzer IRAC and MIPS data combined 
with Herschel PACS and SPIRE data for a uniform sample 
of evolved H II regions shows that the emission in the direc- 
tion of these sources in the range 24 /zm - 500 /xm is well 
represented (for all but three sources) by a 2-temperature dust 
component model with =2. Of these two components, 
the cold one is in the range 20 - 30 K and the warm one in 
the range 50 - 90 K. For one of the sources, G30. 1-0.2, a 1- 
temperature component model with a varying spectral emis- 
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sivity index appears to be a better representation of the ob- 
served SED, although in this case the jijusi (0.32) retrieved 
from the fit is outside the range of commonly accepted val- 
ues. For two other sources, G30.3-0.0 and G30.6-0.1, neither 
a modified blackbody with pdust kept as a free parameter, nor a 
2-temperature component model with fixed (idust agrees with 
the measurements in a satisfactorily fashion. We emphasize 
that, for all the sample sources, the temperatures obtained for 
the warm component are not well constrained due to a lack 
of measurements between 24 /xm and 70 fim. Despite this 
Umitation, these results highlight the importance of includ- 
ing bands shorter than 100 /im for studying sources, such as 
H II regions, which are characterized by multiple temperature 
components. In fact, when the shorter wavelengths are not 
considered in the fitting process, the cold component, as well 
as the warm one, might not be correctly constrained. 

Our study strongly suggests that dust is present in the inte- 
rior of H II regions. Three complementary pieces of evidence 
are found in support of this conclusion: 

1. a statistically significant flux detection, after back- 
ground subtraction, for 8 /im < A < 500 /xm (with an 
average S/N greater or of the order of 3); 

2. a slope for the n^ - d relation slightly steeper than -1.5; 

3. values of IRE > 1 for the entire sample H II regions; 

We emphasize that although the data seemd to favor the ex- 
istence of dust grains in the interior of evolved H II regions, 
the amount of dust in these environments is likely not signif- 
icant. We also remind the reader that, based on our analysis, 
we cannot completely rule out the possibility that the detected 
IR emission is due to background/foreground radiation along 
the Une of sight. 

We argue that partial dust depletion within H II regions ap- 
pears to be caused primarily by radiation-pressure drift, as re- 
cently proposed by Draine et al. (201 1). In this light, we spec- 
ulate that the observed 24 /im emission cannot be attributed 
to a dust population (either VSGs or BGs) co-eval with re- 
spect to the BGs and PAHs located in the surrounding PDR, 
but rather to a new generation of grains, in agreement with 
what was recently suggested by Everett & Churchwell (2010) 
for the case of WBB. 

Finally, we show that far-IR colors can be reliably used for 
identifying unresolved Galactic or extragalactic H II regions 
with respect to younger populations of sources, such as YSOs. 
On the contrary, distinguishing unresolved H II regions from 
SNRs based solely on mid-lR colors appears to be more chal- 
lenging, and may require the additional use of ancillary data, 
such as radio continnuum or RRL. 

The authors would like to thank the anonymous referee for 
his/her comments which allowed us to greatly improve the 
content and presentation of the paper. 
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Fig. 7. — Color-color diagrams obtained by combining PACS (70 ^m and 160 fim) and SPIRE (250 fim, 350 fim and 500 ^m) wavelengths, as well as IRAC 
8 /im and MIPS 24 /im. Shown in the plots are: the 16 evolved H II regions -I- PDRs of this work (blue filled circles), the Hi-GAL point sources from Elia et al. 
(2010) (black filled circles) and the SNRs from Pinheiro-Gon^alves et al. (2010) (black triangles). 
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TABLE 1 
Sample of evolved H II regions 



Name 


1 


b 


^6cm 


R 


D 


^/ 


Te 












"e,shell 


(deg) 


(deg) 


(Jy) 


(kpc) 


(kpc) 


(kpc) 


(K) 


(pc cm"*) 


(pc cm"*) 


(pc cm"*) 


(cm"3) 


(cm"3) 


(cm"3) 




G28.8-0.2 


28.823 


-0.226 


1.47 


4.5 


5.5 




5934 


8.2e4 


4.4e4 




152.7 


95.9 




G29.0-().6 


28.983 


-0.603 


1.01 


5.8 


3.4 




6417 


5.3e4 


2.6e4 


7.0e6 


152.7 


91.3 


5.9e3 


G29. 1-0.0 


29.136 


-0.042 


1.85 


5.8 


11.5 




6417 


2.7e4 


IMA 




43.9 


91.7 




029.1+0.4 


29.139 


0.431 


2.68 


7.4 


13.6 




7012 


3.3e4 


2.3e4 


4.8e6 


41.8 


32.0 


1.7e3 


G29.2-0.0 


29.205 


-0.047 


2.22 


5.4 


10.9 




6268 


3.3e4 


1.6e4 




49.2 


29.1 




G30. 1-0.2 


30.069 


-0.160 


0.95 


4.4 


8.5 




5897 


2.8e4 


2.9e4 




61.9 


62.1 




G30.2-0.1 


30.227 


-0.145 


2.59 


4.4 


6.3* 


8.4* 


9100 


1.0e5 


4.8e4 




147.2/127.5 


83.7/72.5 




G30.3-0.0 


30.277 


-0.020 


0.73 


4.4 


6.2 




5896 


2.5e4 


2.6e4 


2.2e6 


70.5 


72.5 


2.0e3 


G30.4-0.2 


30.404 


-0.238 


3.59 


4.3 


8.0 




5859 


7.0e4 


9.9e4 




90.5 


117.6 




G30.5-0.3 


30.502 


-0.290 


2.96 


4.3 


7.3 




5860 


4.3e4 


7.9e4 




68.7 


108.2 




G30.5-H0.4 


30.467 


0.429 


2.13 


5.7 


3.6* 


11.0* 


6194 


2.0e4 


2.1e4 




59.8/34.2 


61/34.8 




G30.6-0.1 


30.602 


-0.106 


2.40 


4.4 


7.3 




6800 


1.2e5 


8.7e4 




158.2 


123.5 ... 




G30.7-0.3 


30.694 


-0.261 


3.27 


4.5 


8.3 




6800 


1.8e5 


6.7e4 




185.4 


87.3 ... 




G30.8-0.0 


30.776 


-0.029 


62.17 


4.6 


5.7 




7030 


8.6e5 


1.4e5 


4.2e8 


345.2 


88.6 


2.9e4 


G31.0-H0.5 


31.050 


0.480 


1.75 


7.1 


12.9 




6901 


6.3e4 


2.1e4 


5.5e6 


77.8 


34.6 


2.2e3 


G31.1-H0.3 


31.130 


0.284 


1.11 


4.4 


7.3 




5896 


7.4e4 


2.8e4 


7.5e6 


132.7 


64.2 


4.2e3 



Note. — Quoted fluxes are from single-dish measurements by Altenhoff et al. (1979), Downes et al. (1980) and Kuchar & Clark (1997). Galactocentric 
and solar distances are from AB09, based on Lockman (1989) RRL observations. Solar distances denoted with * are computed from radial velocities 
in Lockman (1989) combined with the rotation curve by McClure-Griffiths & Dickey (2007). denotes the far distance solution for sources with 
unresolved distance ambiguity. Details on the derivation of Te, EM and n^ can be found in Section. 3. EM^cm and EM^ocm are, respectively, the emission 
measure computed from single-dish and MAGPIS fluxes, n^^scm and n^ ^Ocm are the corresponding electron densities. EMgheii and n^^jft,,;/ denote the 
emission measure and electron density estimated by assuming that the source is characterized by a shell-like geometry. 
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TABLE 2 
Angular and linear diameters 



Name 


$ 


d 












dshell 




(arcmin) 


(pc) 


(arcmin) 


(arcmin) 


(pc) 


(pc) 


(arcsec) 


(pc) 




2 2 


J.J 


j.\j 


j,\j 




± R 








2.3 


2.3 


3.5 


3.0 


3.4 


3.0 


5 


1 


G29 1-0 


4.3 


14 4 


3.3 


2.1 


11.0 


7.0 






G29 1+0 4 


4.8 


19 


9.7 


3.4 


38.4 


13.4 


12 


0.8 


G29. 2-0.0 


4.3 


13.6 


9.1 


4.1 


28.8 


13.0 






G30.1-0.2 


3.0 


7.4 


3.9 


2.3 


9.6 


5.7 






G30.2-0.1 


2.6 


4.7/6.3 


4.1 


3.5 


7.5/10.0 


6.4/8.5 






G30.3-0.0 


2.8 


5.0 


3.6 


2.1 


6.5 


3.8 


9 


0.3 


G30.4-0.2 


3.7 


8.6 


4.6 


2.1 


10.7 


4.9 






G30.5-0.3 


4.3 


9.1 


5.6 


1.8 


11.9 


3.8 






G30.5+0.4 


5.4 


5.6/17.3 


5.6 


5.1 


5.9/17.9 


5.3/16.3 






G30.6-0.1 


2.3 


4.9 


3.2 


2.3 


6.8 


4.9 






G30.7-0.3 


2.2 


5.3 


4.4 


3.0 


10.6 


7.2 






G30.8-0.0 


4.4 


7.2 


13.8 


8.6 


22.9 


14.2 


6 


0.2 


G3 1.0+0.5 


2.8 


10.5 


5.5 


4.2 


20.6 


15.8 


9 


0.6 


G31. 1+0.3 


2.0 


4.2 


3.1 


3.4 


6.6 


7.2 


6 


0.2 



Note. — The angular diameter 6 and corresponding linear diameter, d, are obtained from the 
single-dish measurements by Altenhoff et al. (1979), Downes et al. (1980) and Kuchar & Clark 
(1997), while Omaj, Smin< ^maj and Amin are derived from the MAGPIS data (Section 3 and Figure 7). 
All linear diameters are estimated with the heliocentric distances quoted in Table 1. 6 shell and A^heii 
denote the shell angular size and linear thickness in the hypothesis of a shell-like geometry of the 
source. 
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TABLE 5 
Best-fit SED parameters 



Name 


dust, 2 




Pdust 






CDFi, 




(K) 


Tw,3 


T„3 


X2(^ 


CDFj, 




(K) 


(K) 








(K) 


(K) 


(K) 








72.66 


10.26 


1.2x 10"^ 


0.33 


9.58 


0.98 


17.73 


0.50 


55.92 


1.25 


0.86 


0.35 


ozy.u-u.o 


Q1 
DD.iSl 


y.o / 


J. OX lU 


U.3D 


o.y3 


u.y / 


21. Z/ 


U.iSU 


D4.1 / 


u.yo 




\J.DO 


029.1-0.0 


58.98 


13.89 


7.3x 10~' 


0.73 


16.40 


1.0 


29.46 


0.49 


96.31 


9.52 


1.54 


0.54 


Cj29.1-t-0.4 


60.96 


7.65 


1.8x 10~^ 


0.41 


22.93 


1.0 


23.40 


0.90 


57.12 


15.59 


0.35 


0.16 


029.2-0.0 


57.15 


4.82 


3.7 X 10~^ 


0.32 


24.29 


1.0 


25.62 


0.40 


98.06 


1 1.83 


0.41 


0.19 


G30. 1-0.2 


65.04 


2.14 


3.2x 10"' 


0.09 


2.91 


0.59 


19.81 


0.89 


51.49 


0.72 


2.34 


0.69 


G30.2-0.1 


82.45 


1.12 


1.0x10-2 


0.11 


7.34 


0.94 


20.63 


0.89 


60.14 


1.81 


1.11 


0.42 


G30.3-0.0 


55.01 


8.79 


1.1 


0.50 


24.49 


1.0 


31.91 


1.01 


98.37 


10.54 


7.16 


0.97 


G30.4-0.2 


65.56 


16.24 


4.6x ur' 


0.67 


17.76 


1.0 


26.70 


0.93 


70.55 


12.00 


0.15 


0.08 


G30.5-0.3 


56.93 


12.37 


4.2x ur' 


0.62 


11.83 


0.99 


24.65 


1.10 


53.00 


11.69 


0.92 


0.37 


G30.5+0.4 


57.25 


7.00 


8.0xHH 


0.38 


8.73 


0.97 


18.95 


0.55 


48.55 


0.98 


0.40 


0.18 


G30.6-0.1 


22.60 


1.36 


2.75 


0.21 


102.20 


1.0 


27.95 


0.38 


95.07 


11.84 


0.79 


0.98 


G30.7-0.3 


56.15 


1.83 


1.31 


0.11 


9.11 


0.97 


36.15 


0.92 


97.33 


10.79 


2.98 


0.77 


G30.8-0.0 


51.61 


4.63 


8.2x10-' 


0.25 


10.29 


0.98 


25.94 


1.08 


50.33 


16.81 


0.40 


0.18 


G3 1.0+0.5 


56.84 


3.58 


6.0 X 10-2 


0.28 


28.07 


1.0 


22.58 


0.60 


56.10 


7.18 


1.22 


0.46 


G31. 1+0.3 


68.16 


18.55 


3.0X 10-5 


0.71 


4.28 


0.77 


19.72 


0.64 


52.33 


1.02 


0.23 


0.11 



Note. — Results obtained by fitting the 24 /xm < A < 500 /xm data points. Tdust,2< '^T^^sti Xii^ are the best-fit 
parameters and corresponding for a 1-temperature component model with variable ^dust- Tc,3, cr^j, T„_3, aT„, and X2»^ 
are the best-fit parameters and f™ a 2-temperature (warm and cold) component model with /3^„s, equal to 2. PDFi, and 
PDF2, are the probability density functions associated to each model. 



TABLE 6 
Luminosities and IR excess 



Name 


(s-') 


log(L(Ly a)/L©) 


log(Lre,i/7//L0) 


log(L;/e_PD«+H///L0) 


JREhii 


JSEpoR+Hii 


G28.8-0.2 


48.67 


4.32 


4.37 


5.22 


1.01 


1.21 


G29.0-0.6 


48.07 


3.7 


4.5 


4.95 


1.21 


1.34 


G29.1-0.0 


49.39 


5.05 


5.22 


5.74 


1.03 


1.14 


G29. 1+0.4 


49.69 


5.34 


5.55 


6.14 


1.04 


1.15 


G29.2-0.0 


49.43 


5.08 


5.64 


5.97 


1.11 


1.17 


G30.1-0.2 


48.86 


4.51 


4.71 


5.13 


1.04 


1.14 


G30.2-0.1 


48.95/49.20 


4.60/4.85 


3.46/3.71 


5.29/5.54 


0.75/0.76 


1.15/1.14 


G30.3-0.0 


48.47 


4.13 


4.52 


4.92 


1.09 


1.19 


G30.4-0.2 


49.38 


5.04 


5.96 


6.29 


1.18 


1.25 


G30.5-0.3 


49.22 


4.88 


5.21 


5.66 


1.07 


1.16 


G30.5+0.4 


48.46/49.43 


4.11/5.08 


4.61/5.57 


4.73/5.69 


1.21/1.10 


1.15/1.12 


G30.6-0.1 


49.10 


4.75 


5.09 


5.36 


1.07 


1.13 


G30.7-0.3 


49.35 


5.00 


5.55 


5.96 


1.11 


1.19 


G30.8-0.0 


50.29 


5.95 


6.62 


6.89 


1.11 


1.16 


G3 1.0+0.5 


49.46 


5.11 


5.86 


6.01 


1.15 


1.17 


G31. 1+0.3 


48.79 


4.45 


4.80 


5.23 


1.08 


1.17 



Note. — Lyman continuum flux (Nl^c). Lya Imninosity (L(Ly a)), IR limiinosity (Lir^hii and 'Lir^pdr+hii) and IR 
excess (JSBhii and ^SBpbr+hu) for the 16 H II regions (and associated PDRs) in our sample. 
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Fig. A1. — 3-color images of the sample evolved H II regions. Left panel: MAGPIS 20 cm (red), MIPS 24 /im (green), IRAC 8 fim (blue) data. Dashed 
lines indicate cuts in latitude and longitude operated to create the profiles in Appendix B. Right panel: MAGPIS 20 cm (red), SPIRE 250 fim (green), PACS 70 
fim (blue) data. G29. 1-0.0 and G29. 2-0.0 are shown in the same panel (third row). Areas enclosed by white dashed lines denote approximate apertures using 
MAGPIS 20 cm data. 
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Fig. B1. — Longitude (right) and latitude (left) profiles, obtained by slicing MAGPIS 20 cm, IRAC 8 fim, MIPS 24 fim and SPIRE 250 /im postage stamp 
images of the sources. All images are convolved to SPIRE 250 fim resolution (18") prior to slicing. 
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